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The temperature-dependent nonlinear conductance for transport of a Luttinger liquid through a barrier is calculated in the 
nonperturbative regime for g — l/2 — e, where g is the dimensionless interaction constant. To describe the low-energy behavior, 
we perform a leading-log summation of all diagrams contributing to the conductance which is valid for |e| <C 1. With increasing 
external voltage, the asymptotic low-temperature behavior displays a turnover from the J^/s— 2 to a universal T 2 law. 

PACS numbers: 72.10.-d, 73.40.Gk 
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I. INTRODUCTION 

The interplay between electron-electron interactions 
and disorder has attracted a great deal of recent in- 
terest in condensed-matter theory. Many-body correla- 
tions play an essential role in quasi-one-dimensional (ID) 
fermionic systems, where the usual Fermi liquid behavior 
is destroyed by the interactions. The generic behavior of 
many interacting ID fermion systems isjpstead described 
in terms of the Luttinger liquid modehlllu One may then 
ponder how transport in such a ID "quantum wire" is 
affected by the presence of impurities or tunnel barriers. 
A quantitative answer to this question is of immediate in- 
terest for several experimental setups, e.g. the tunneling 
of edge state excitations in the fractional quantum Hall 
(FQH) regimejj'cl or transport in narrow high-mobility 
heterostructure channels BEl 

The effect of one or a few barriers on the conductance 
of a Luttinger liquid has first been studied by Kane and 
Fisher (KF).Q The Luttinger liquid model captures the 
low-energy properties of correlated electron transport if 
two conditions are met: one has short-ranged electron- 
electron interactions (i.e. the long-range 1/r tail of the 
Coulomb potential is screened by nearby gates), and 
electron-electron backscattering (BS) merely constitutes 
an irrelevant perturbation.@ For the spinless case, BS is 
the exchange process of forward scattering, and therefore 
it is included readily by a redefinition of the interaction 
constant. All effects of the electron-electron interaction 
arc then encompassed in a single dimensionless parameter 
g, which is approximately given by g » (1 + U /2Ep)~ 1 / 2 . 
Here, U measures the interaction strength and Ep is the 
Fermi energy. We will consider the case of repulsive inter- 
actions, where one has g < 1. The noninteracting Fermi 
liquid case is recovered for g = 1. An ideal realization of 
this model is found in the tunneling of edge state excita- 
tions in the FQH regime. As shown by Wenj3 these exci- 
tations are described by a (chiral) Luttinger liquid with 
g = v, where v is the filling factor. Since v is a topo- 



logical quantity controlled by the bulk properties of the 
FQH bar, this provides an excellent experimental testing 
groundo In addition, recent transport experiments in ID 
quantum wiresa have revealed Luttinger liquid behaviors 
directly for correlated electrons. 

It has been demonstrated that the presence of an im- 
purity potential term (the simplest realization of disor- 
der) leads to a metal-isolator quantum phase transition at 
T = when going from attractive to repulsive electron- 
electron interactions .13 For g < 1, any (whatsoever weak) 
barrier will therefore result in zero linear conductance, 
G = 0. One is then concerned with the low-energy prop- 
erties of the conductance, i.e. the behavior of G(T, V) in 
the limit T, eV/ks <C 2k, where the Kondo temperature 
Tk sets the appropriate scale (see below). As indicated 
by the RG flow, this is a nonperturbative problem for 
weak barriers, and the past few years have shawn_cpnsid- 
crablc effort being devoted to this problem. BOtTlI 2 ] 

By using a simple scaling argument, Kane and Fisher 
argued that the linear (V — > 0) conductance should 
behave ,-as T 2 /f~ 2 in the asymptotic low-temperature 
regime.Q They obtained this law from matching the RG 
flow in the nonperturbative regime onto the large-barrier 
perturbative flow, thereby exploiting universality argu- 
ments. On the other hand, any finite voltage leads to 
a ^-independent (universal) T 2 power law instead of the 
T 2/ g -2 behavior. As one can never be sure about the va- 
lidity of scaling arguments, we believe it is of use to have 
exact results for the nonequilibrium finite-temperature 
transport as provided here. We note that the exact solu- 
tion for the case g = 1/2 lias already been given within 
the bosonized description£3 The linear conductance for 
g = 1/2 has also been calculated within an equivalent 
fermionic model by Kane and Fisher ,0 who extended an 
earlier zero-temperature approach by GuincalU to finite 
temperatures. Since the laws T 2 / g ~ 2 and T 2 coincide 
for g — 1/2, however, it is crucial to consider g ^ 1/2. 
In this paper, we start out from a dynamical approach 
and present a nonperturbative leading-log summation for 
G(T,V) in the range g = 1/2 - e with |e| < 1. We 
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systematically calculate the conductance in form of the 



G(T, V) = a 3 e j ln J [max(fc B T, eV)] 



(1.1) 



Our study is related to recent work by Matveev, Yue 
and Glazman.tJ They have carried out a leading-log sum- 
mation for weak electron-electron interactions. For the 
Luttinger liquid, this would correspond to the case g = 
1 — e with |e| <C 1. Since the Kondo temperature van- 
ishes with an essential singularity as g approaches 1 from 
below, it is desirable to have similar nonpcrturbative re- 
sults for stronger electron-electron interactions. Their 
technique relies crucially on the assumption of weak in- 
teractions and can therefore not be employed directly in 
other regions of parameter space. 

Basically parallel to our work, a powerful approach 
based on the thermodynamic Bethe ansatz has been put 
forward by Fendley, Ludwig and Saleurtil who studied the 
linear conductance at g — 1/3. They mapped the orig- 
inal boson model of the Luttinger liquid .onto a mass- 
less odd boson model which is integrable.Ej A basis of 
massless fermion-type quasiparticles can be derived, in 
which the S matrix for impurity scattering can be calcu- 
lated exactly. The distribution function and the density 
of states of these quasiparticles are given in terms of inte- 
gral equations with kernels determined by the completely 
elastic and factorizable scattering matrix of the bulk. Ex- 
pressed in terms of these quantities, a Boltzmann-type 
rate expression for the conductance may then be derived. 
Switching to a suitable momentum integral representa- 
tion, our solution for the current obtained from a dy- 
namical approach can be transformed such that we can 
read off the S matrix for impurity scattering, the pseu- 
doenergy and the density of states of the quasiparticles 
for g = 1/2 — e and finite voltage. Therefore, our dynam- 
ical approach provides an independent check for results 
obtained by the thermodynamic Bethe ansatz, and fur- 
thermore gives a nontrivial prediction for the density of 
states. 

The outline of this paper is as follows. In Sec. II, we 
present the model for our study and give the formal so- 
lution. The general diagrammatic expansion is explained 
in Sec. Ill, and the leading-log summation based on this 
diagrammatic expansion is discussed in great detail in 
Sec. IV. Results for the asymptotic low-energy behavior 
as well as a conjecture for the exact solution are presented 
in Sec. V, followed by concluding remarks in Sec. VI. 
Some technical details of the leading-log calculations are 
given in the Appendix. 



II. GENERAL FORMALISM 

To describe the problem of correlated electron tr; 
port, we employ the standard bosonization method.EJ 



The fermion field operator is expressed in terms of two 
bosonic fields <fi(x) and 9(x), 

ip\x) ~ exp[m(v / 7r^(a;) + kpx) + iy/TT(j)(x)] , 

n odd 

where the boson fields obey the commutation relation 

[cl>(x),6(x')} = -(i/2)8gn(aj - x') . 

Taking short-range electron-electron interactions and ne- 
glecting backscattering, this leads to the generic Lut- 
tinger liquid Hamiltonian (we put % = 1) 



H 



dxlgid^) 2 +g-\d x ey 



(2.1) 



where g is the dimensionless interaction constant and v 
denotes the sound velocity. The Fermi velocity is then 
given by vp = gv. For the physical case of repulsive inter- 
actions, one has g < 1. Although this model is universal, 
actual computations, e.g. of correlation functions, neces- 
sitate introduction of a non-universal cutoff parameter w c 
which is related to the bandwidth of the non-interacting 
problem. In the end, -since Hq is quadratic, everything 
can be solved exactly.El n 

To model the impurity, we follow KF0 and consider 
a short-ranged scattering potential centered at x = 0. 
Omitting irrelevant multi-electron backscattering pro- 
cesses, the important contributions are captured by the 
2fcF-component Vq of this scattering potential, and one 
has 

H = H + V cos[2^(9(0)] + eV (9(0) /y/n . (2.2) 

We have also included an external voltage V by assuming 
the voltage drop to occur at the barrier. This is well 
justified for no t to o small barriers. 

The model (2.2) describes scattering by an impurity 
potential with effective strength Vo. Interestingly, results 
for a very high .barrier Vq S> w c can be obtained via 
an exact dualitycj from the case Vq <C u> c . For g < 1, 
the large-barrier problem can be treated in lowest order 
in the number of tunneling transitions for any T and 
V, whereas in the weak-barrier case perturbation theory 
in Vq is not possible at sufficiently low temperature and 
small voltage. Our subsequent treatment deals with the 
interesting nonperturbative regime. Putting g = 1/2 — e, 
our results presented in Sec. V hold under the conditions 

Vq/u c « 1 , |e| « 1 , eV/k B T K « 1 , T/T K « 1 , 

(2.3) 



where Tk is the Kondo temperature 

t k = (MK) s/(1 ^ 9) M/t B 



(2.4) 



taken at g = 1/2 — e. The Kondo temperature provides 
an estimate for the extent of the nonpcrturbative regime. 
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Note that 7k — > for g — > l - as alluded to in the Intro- 
duction. 

In the following, we study the nonlinear static conduc- 
tance G(V, T) — dl/dV. From the bosonized form of the 
current operator, it takes the form 



G(V,T) = (e/v^F) (d/dV) \im (0(0, i)) , 

t — >oo 



(2.5) 



wh ere (• • •) denotes the thermal average over all mod es 
of ( |2.2| ) away from x = 0. Remarkably, the model (2.2) is 
identical to the prohlmi xrf.a Brownian particle moving in 
a periodic potentialEj'oESThis can be directly shown by 
a unitary transformation,!!-!! and we exploit such a map- 
ping here. The (9-field taken at the location of the im- 
purity, #(0), corresponds to the position operator of the 
Brownian particle, and the Ohmic dissipation is provided 
by the plasmon modes away from the barrier. Thereby 
the conductance can be expressed in terms of the nonlin- 
ear mobility of the Brownian particle which is coupled to 
an Ohmic heat bath. 

An exact formal expression for the conductance can be 
obtained from a real-time path-integral approach due to 
Feynman and Vern on £3 Employing the unitary transfor- 
mation of Eq. (2.2) onto the tight-binding limit of Brow- 
nian motion in a periodic potential, one can write this ex- 
pression in a "Coulomb gas" representation of in tera cting 
discrete charges.Q Alternatively, our expression ( [2.6] ) with 
(2.7) can be found by expanding the impurity propaga- 
tor in terms of auxiliary charge pathsJl3 The 6(0) paths 
can then be eliminated analytically and one has to sum 
over Hubbard-Stratonovich paths instead. These auxil- 
iary paths can be parametrized by the discrete charges 

In the end, the exact formal expression for the conduc- 
tance takes the form of a power series in Vq , 



G(V, T)/G = 1 - (27r/ e y) Im U(V, T) , 



(2.6) 



where Go = ge 2 jh and U describes the interacting charge 
gas 

oo °° 

U(V, T) = £ (iV ) 2m / d Tl dr 2 ■ ■ ■ dr 2m ^ (2.7) 
m=l 





2m-l 



with the interaction factor 



2m 

IF,,, = exp | t,J S ( T jk)£,k 
\j>k=l 



(2. 



The interaction potential between the charges is given by 

S(t) = 2gln[{uj c /Trk B T) sinh(7rfc B Tr)] . (2.9) 

The 2m — 1 integration times Tj in order Vg m are the in- 
termediate times between the 2m successive charges, and 



Tjk gives the time interval between charges and 
Charge configurations {£} contributing to (2/7) are re- 
stricted as follows. (1) In every order m, the 2m charges 
have zero total charge, Ylj=i £j = (neutrality condi- 
tion). (2) The cumulative charge quantities Pj t m 



Pj, 



2 m 

i=J+l 



■Eft 



can never be zero since otherwise the corresponding phase 
factor sin(irpj jm g) would vanish. Expressed in graphical 
terms, there are no diagrams that can be subdivided into 
disjoint parts each of them representing a neutral config- 
uration. (3) By convention, the first charge is = — 1. 
With that, all Pj, m are positive integers. 



(a) 













(b) 



■ P = 2 



■ P=l 



.p = 



FIG. 1. Charge configura tion ( h+) contributing in 

order m = 2 to the expression (2.7). (a) Diagram representing 
this configuration. Time flows along the horizontal line, ver- 
tical lines symbolize charges £ = ±1. (b) Accumulated charge 
Pi = X/j>i £i during the course of the path. Circles stand for 
the momentary position of the path. For the shown config- 
urati on, a n overall phase factor 2ne is acquired according to 
Eq. fl210|). 



To illustrate these rules and our subsequent diagram- 
matic approach, Fig. [I] shows a charge c onfig uration for 
m = 2 which does contribute to Eq. (2/7). Put ting 
g = 1/2 — e, the phase factors appearing in Eq. (2.7) 
can be written as 



sin(Trpg) 



(_l)(p-l)/2 , podd, 

irep (— 1)(p -2 )/ 2 , p even. 



(2.10) 



As one moves away from p — 0, one has to pay the 
phase factors 1, 2-7re, —1, — Aire, . . . for dwelling on states 
labelled _p — 1, 2, 3, 4, . . ., respectively. This is displayed 
in Fig. |l](b). The charge configuration in Fig. [y has a 
phase prefactor 27re. On the other hand, the breathing 
mode integral over the interior dipole [the T2-integration 
of the (m = 2)-term in Eq. ([2.7j)] has a 1/e singularity 
due to the r~ 1+2c short-time behavior of the intradipole 
interaction exp[— S(t)]. Hence the combined expression 
is well defined even when e = 0. 
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The essential scale which separates the nonperturba- 
tive from the perturbative regime may ba-derived by us- 
ing an adiabatic renormalization scheme.ta Rather than 
working in frequency space, we take advantage of a time 
formulation. To obtain a renormalized frequency scale 7, 
we integrate the length of a dipole from zero up to 1/7, 
where the inverse time scale 7 is self-consistently given 
by the short-time part of the breathing mode integral. 
We have 



/•1/7 

7 = 2 7 reV 2 / dr 
Jo 

From this we obtain 

7 = 7 (w c /7) 2e 



1 



irVn I LOr 



Kr) 1 - 



2 c 



(2.11) 



where 7 = ttVq /oj c 



The frequency 7 is the effective inverse time scale of the 
problem, and 7 is the corresponding quantity when e = 0. 
The related temperature 7 //cb coincides with the previ- 
ously introduced Kondo temperature (2.4) for g = 1/2— e. 

To develop a general computational strategy, it is con- 
venient to split up every time integral over the length of 
a dipole into a contribution from the short-time regime 
< t < I/7 and a residual contribution from lengths 
t > I/7. We shall refer to the short-time part as col- 
lapsed dipole. The crucial observation guiding us towards 
this terminology is as follows. Within leading-log ac- 
curacy, a collapsed dipole has zero dipole moment and 
therefore no interactions with other charges. This prop- 
erty allows us to carry out the grand canonical sum over 
all possible arrangements of collapsed dipoles between 
two confining charges in an exact manner. The strategy 
is then completed by developing a systematic scheme to 
calculate the contributions from all the residual time in- 
tervals Tj > I/7 in the asymptotic low-energy regime 



(2.3) 



There are two types of collapsed dipoles, namely (H — ) 
and (— +). If they describe hops forth and back from 
states labelled by an even p value, the contributions of 
these two dipoles cancel each other completely. This is 
caused by the different sign of the phase factors — irep 
and Trep which are connected with the dipoles ( — h) and 
(-1 — ), respectively. Hence there is no insertion of a gas of 
colla psed dipoles for the even time intervals T2,T4, . . . in 
Eq. (^7j), where the path is in an even p state. However, 
if the dipoles describe hops forth and back from states 
labelled by odd p, the dipoles (— +) and (H — ) come with 
phase factors 7re(l +p) and 7re(l —p), respectively. Note 
that the second type does not contribute when p = 1. 
The sum of both contributions is independent of p. Al- 
lowing next that the collapsed dipole ( — h) [as well as 
(H — )] moves freely within an interval r, we find a to- 
tal contribution — jt. Finally, if we consider instead of 
the two dipoles a grand canonical noninteracting gas of 
these being confined within an odd-time interval t, the 
total contribution is summed to an exponential factor 
exp(— tt). 



To summarize these findings, the overall effect of the 
collapsed dipoles is to cause exponential suppression of 
every long odd time interval on a length 1 /% while long 
even time intervals are not suppressed on this scale. In 
the diagrammatic expansion given below, we indicate the 
insertion of a grand canonical gas of collapsed dipoles by 
a square box. We shall refer to bare diagrams when they 
are free of collapsed dipoles and to dressed diagrams when 
every odd time interval is dressed by a square box. 



III. DIAGRAMMATIC EXPANSION 



In the following 
sion (2J3) with (|2 



we wish to evaluate the exact expres- 
) for the nonline ar co nductance in the 
low-energy regime specified in Eq. (2.3). In that 
leading-log summation becomes possible around g = 1/2. 
With g — 1/2 — e, we compute the coefficients dj in the 
series 



00 



G(V, T) = ^ dj (V, T) e } hr" [max(fc B T, eV)] 

3=0 



(3.1) 



Contributions ~ e In 5 [max(/cBT, eV)] with j < k are dis- 
regarded in our calculation. This approximation is valid 
under conditions ( |2.3| ). In this section, computation of 
the functions dj (V, T) is presented in detail. 



FIG. 2. Diagram giving the contribution Ui to the con- 
ductance. The square stands for a factor exp[— 77-] due to a 
gas of collapsed dipoles. 

We start by rewriting the quantity U ( V, T) defined in 
( |2.7| ) in terms of a diagrammatic expansion, 



U(V,T) = X>«- 



(3.2) 



71=1 



The quantity U n is the sum of all arrangements of 2n 
charges complying with the rules stated above. The pro- 
cedure is then as follows. (1) We have to draw all possible 
bare diagrams with In charges of total charge zero which 
cannot be subdivided into disjoint neutral parts, and the 
first charge has to be £1 = — 1. (2) In the next step, the 
bare diagrams are dressed by decorating all odd time in- 
tervals T2j-i (j = 1, . . . , n) with square boxes. Each of 
the square boxes symbolizes insertion of a grand canoni- 
cal gas of collapsed dipoles, i.e., a factor of exp[— 775^-1]. 
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(I) 



FIG. 3. Diagram giving the contribution Ui- 
As outlined before, there are no such insertions for even 

(II) 

time intervals T2j . To avoid double counting, integrations 

over even times are free of collapsed contributions, and 
therefore we restrict these integrations to *fT 2 j > 1. The 
term U n contributes to order e™ -1 (and to higher orders). 
To illustrate this, the diagrams representing the terms U\ 
to U4 are shown in Figs. § to g. While in orders n = 1 
and n = 2 only one diagram is found, respectively, there 
are two diagrams for n = 3 and five diagrams for n = 4. (Ill) 

(!) ~T 



(II) 



(IV) 



FIG. 4. The 2 diagrams yielding the contribution U3 



The dressed diagram U\ yields already the exact solu- 
tion for g = l/2 since all o ther contributions are at least 
of order e. From Eq. (2.7), we find 



Ui = -V 2 / dr exp[-(7 + ieV/2)r - S(t)} , 
Jo 

which gives indeed the exact solution for g = 1/2 in terms 
of the digamma function if>(z), 



(3.3) 



Thus the notion of collapsed dipoles provides a-xemark- 
ably simple derivation of this important result .113 For fi- 
nite e, the quantity U\ is the same as for g = 1/2 within 
leading-log accuracy, except for 7 being replaced by 7. 
At this point, we would like to remark that the dia- 
gram U\ gives the correct limiting behavior G/Gq — ► 
as T — > 0, V — * for any Vq. Hence all higher-order 
dressed diagrams should give vanishing contributions in 
the zero-energy limit, and this indeed we find. 



(V) 





































FIG. 5. The 5 diagrams contributing to U4. Diagrams (I) 
- (IV) are of type (A), and diagram (V) is of type (B). 

Generally, besides the lowest-order diagram in Fig. ^ 
we may distinguish two types of dressed diagrams. Dia- 
grams of type (A) are of the form ( X + +). Here, 

X stands for insertion of any arrangement of extended 
dipoles (H — ) or £]-+) between the outer double charges 

( ) and (++)£j Since there are two types of extended 

dipoles, we have 2 n ~ 2 dressed diagrams of type (A) con- 
tributing to order U n . These dipoles are of finite length 
in contrast to the collapsed dipoles considered before. 
Hence they are interacting with each other and also with 
the outer charge pairs. However, because of the insertions 
of collapsed dipoles (indicated by the squares in the dia- 
grams) , they are narrow compared to typical inter-dipole 
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distances. All other diagrams contain at least four pairs 
of equal charge. We shall refer to them as diagrams of 
type (B) . The simplest one of this type contributes to U4 
and is diagram (V) in Fig. @. Remarkably, up to order 
U3 only diagrams of type (A) are encountered. 

Now we have to extract and evaluate the leading loga- 
rithmic contributions in the terms with n > 1. It simpli- 
fies notation to use scaled times and energies henceforth, 



= 27rfc B T/7 , v = 2igeV/^ . (3.4) 



Xj = ~fTj , 



According to the conditions ( |2.3| ), we are interested in 
the asymptotic regime 



u < 1 



M < 1 



(3.5) 



Now every odd time integration J dx 2 j +1 has a weight 
function exp(— x 2j+1 ) due to the insertion of collapsed 
dipoles while all even time integrations are free of this 
weight. Hence it follows that the expression U n is domi- 
nated by the regime 



X 2j+1 ^ X 2k 



J' 



0,1,. 



1 



(3.6) 



and the logarithms we wish to extract must arise from the 
integrations over the long even times. In the following, we 
describe in detail the general strategy we have employed 
to evaluate diagrams like the ones shown in F igs. || 0. 

We first rewrite the interaction factor (2.8) using the 
definitions 



2n 

M n = exp ^2 ijR{xjk)ik , 

j>k=l 

2n£ M" 2e , 



N n = (2/u) 
R(x) — sinh(ux/2) 



such that 



W n = (7K) (1 - 26)n (u/2) n M n N n 



(3.7) 

(3.8) 
(3.9) 



(3.10) 



The factor M n includes the interactions of 2n charges for 
g = 1/2, while the factor N n encapsulates the remaining 
e-dependent part of the interactions. The usefulness of 
the split-up ( 533 ) is demonstrated now. 

Having drawn a certain diagram contributing to U n , 
say Un , we take advantage of a general deco mposition 
theorem for the interaction factor M n in (3.7). The de=. 
composition theorem corresponds to Wick's theoremEj 
for fermionic fields and originates from the equivalence 
of the interaction factor exuL— S(t)) for g = 1/2 with a 
free ID fermion propagator J13 By virtue of this theorem, 
we can decompose M n — which contains n(2n — 1) pair 
interactions — into a sum of n! terms containing only 
n intrapair interactions of neutral pairs. This procedure 
results in the corresponding decomposition o f a gi ven nth 
order diagram into n! graphs. The split-up (3.10) allows 
us therefore to topologically reduce the diagrams to graphs 
which are much easier to evaluate. Thereby the quantity 



U n takes the form of a sum over all diagrams (labelled 
by an index (j)) where each diagram is represented by a 
sum of n! graphs, 



U„ 



,71 — I n; 

e u 7 \^X^TT(j,r) 



2tt 



j r=l 



(3.11) 



where U„ describes the contribution of the rth graph 
to the nth-order diagram j under study (with the pref- 
actor chosen accordingly). Pictorially, one obtains the 
graphs by grouping the 2n charges into n neutral pairs. 
Clearly, there are n! different assignments and hence n! 
different graphs for each diagram, with the sign of each 
graph given by the number of crossings of lines connect- 
ing the paired charges. Fig. ^ shows the 2 graphs for the 
diagram in Fig. || resulting in U%, and in Fig. ^ the 12 
graphs due to the two diagrams of U3 are displayed. 

In the final step, we group all graphs Un'^ into classes 
with respect to their particular dependence on the even 
time intervals x 2 j- Each of these classes is then evalu- 
ated in leading-log accuracy, and the remaining odd time 
integrations can be carried out without further approxi- 
mation. 

So far we have not taken i nto account the residual in- 
teraction factor N n given in (3.8). We have to treat this 



term separately as it does not fit into the decomposition 
scheme. Under condition ( |3.6| ) and within leading-log ac- 
curacy, this factor has the same form for all diagrams of 
type (A) and is 



N n A) = (x 2 + x 4 + ---+x 2n _ 2 ) 



8c 



(3.12) 



Unfortunately, for graphs of type (B), it is not so straight- 
forward to include the factor N n in the explicit calcula- 
tion. Therefore, in this paper, we confine ourselves to 



the regim e where it is consistent to put N n 
Sec. IV C and Appendix). 



(B) 



1 (see 



IV. LEADING-LOG EVALUATION OF THE 
CONDUCTANCE 

Following the general diagrammatic approach outlined 
in the previous section, we have explicitly calculated the 
leading logs in the expressions U n for n < 4. For larger 
n, the explicit calculation of all diagrams becomes pro- 
hibitively cumbersome. However, for the asymptotic low- 
energy corrections only the restricted set of diagrams of 
type (A) contributes, as we shall explain below. This per- 
mits an exact evaluation of the asymptotic corrections. 
The corresponding summ ation of all diagrams of type 
(A) is carried out in Sec. V A . In the remainder of this 
section, we illustrate the computation of U n for n < 4. 
That will also lead us towards the exact solution of this 
problem, see Sec. 



G 



A. Contribution U2 

To proceed according to the strategy of Sec. Ill, we 
have to first draw all diagrams in order n = 2. There is 
only the single diagram shown in Fig. ||, which is of type 
(A) and contributes in order e (and higher orders) to the 
series (3.1). In a second step, this diagram may be topo- 
logically reduced into the two graphs shown in Fig. ||, as 
follows from the decomposition of the interaction factor 

]\/[ 2 — e R(x 1 )+R(x 3 )~R(x 2 )~R{x 1 +x 2 )~R(x 2 +x 3 )-R(x 1 +x 2 +x 3 ) 
_ -Rix^-RiXi+Xri+Xz) _ -i?(x 1 +a; 2 )-_R(2; 2 +X3) 



where R(x) is given by Eq. (Si)). Both graphs fall into 
different classes, i.e. they have a different functional de- 
pendence on the lo ng ev en ti me int erval X2- In particular, 
according to Eqs. (2.7) and (3.11), we have to evaluate 
the terms 



U, 



(i) _ 



dxidx; 



dX2 



x 8e e -[l+v/2](xi+x 3 )~vx 2 



sinh[?i(a;2 + X\ 
dx\dx% 



Xz)/2\ sinh[uiE2/2] ' 
dx 2 

x 8e g— [l+v/2](x 1 +x 3 )-vx 2 
X - 

sinh[u(a;2 + Xi)/2] sinh[u(a;2 + 

where we have used ( |3 . 1 2[ ) . Recalling that the dominant 
contributions to the triple integrals come from the region 
x\ , a; 3 <C X2, it is clear that the logarithms which we wish 
to extract arise from the ^-integrals. To proceed, we use 
the approximation 



it exp[uir + u(a + b)/2] 
4sinh[u(a; + a)/ 2] sinh[u(ir + b)/2] 

(1 _ e — u(6-a)\-l ^ _ e u(b-a)\ 



(4.1) 



x + a 



where x corresponds to the long even time x 2 and a, b are 
linear combinations of short odd times x\, x%. The terms 
neglected on the r.h.s. are subleading as they do not give 
logarithmic contributions. It is at this point where we 
invoke the leading-logarithmic approximation. 

Identifying a and b with q, the contributions to the 
integral over the even time are in the form 



„8c 



dx- 



<1 



-(u+v)x _ e q(u+v) £ 2 ( e ) 



(4.2) 



Since the odd time integrations produce a purely imagi- 
nary contribution, we need to consider only the real part 
of L2(e), i.e. we may replace u + v by \u + v\ in this term 
when performing these integrations. Within leading- log 
accuracy, we have 



L 2 (e) = 



dx- 



„l-8e 



= -lnl 



;| (\ + 0{e\n\u 



(4.3) 



(2) 



r 



7=r 



FIG. 6. The two graphs obtained from the decomposition 
theorem for the contribution U2 shown in Fig. ^. The relative 
signs are + for graph (1) and — for graph (2). The curves 
display interactions between paired charges. 



Upon using Eqs. (4.1) and (4.2), the resulting expres- 
sions for C/ 2 (1) and may be combined in a form which 
will find straightforward generalization to higher orders, 



Imy^L^^ L 2 (e) / dxidx^e 

u Jo 



-(xi+x 3 ) 



(4.4) 



E 

{^■=±1} 



CT103 



sin[|?;|(o-i3;i + 0-3X3)72] 
sinh[w(cria;i + o- 3 a;3)/2] 



The remaining odd time integrations pose no problem 
and are done easily. 

Collecting all terms together, we find the leading-log 
contributions to the conductance resulting from U2 in all 
orders in e, 



lmU 2 = 4e-L 2 (e)Im 

7T 



+trV 



2u 



2 it 



(4.5) 



The function ip' is the derivative of the i/j function. To- 
gether with the contribution due to Ui, this gives all 
terms in order e (plus partial summation of higher-order 
terms in e). 



B. Contribution U3 

To proceed further, let us now study the next order. 
The two diagrams contributing to order t/3 are displayed 



7 



in Fig. ||. Applying the before-mentioned decomposi- 
tion scheme for the interaction factor M3, we find the 
12 graphs shown in Fig. Again we group them into 
classes regarding their dependence on the even intervals 
X2,X4. The contributions U^' (j = 1, 2; r = 1, . . . , 6) 
corresponding to the 12 graphs in Fig. fj] are expressed as 



U. 



(l,r=l,...,6) _ 



dxidx^dx^ 
Ji 



dx2dx± 



x e 



-[l+v/2](xi+X3+x 5 )-vx 3 p -v{x2+Xi) jy(A) j(l,r) 



(2,r=l,...,6) 



dx\dxj,dx^ J dx2dxi 



-[l+v/2](x 1 +x 3 +x 5 ) e -v{x 2 +Xi) jy(A) A2,r) 

3 J 3 



where Nn is given in ( [3.12 ), and where each ff' r ' be- 
longs to one of the following two classes, 

= (smh[u(x 2 + 24 + o)/2] 
x sinh[u(a;2 + x A + b)/2] sinh[wc/2]^ , 
/f 5 = (smh[u{% 2 + o)/2] 

x sinh[u(a;4 + 6)/2] sinh[u(a;2 + x 4 + c)/2]^ 

Here, a, 6, c denote certain linear combinations of the 
short odd time intervals £1,0:3,2:5 <C X2,a;4- This classi- 
fication is again motivated by the fact that logarithmic 
terms are attributed to the domain of long even time 
intervals. 
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where we have reintroduced = x 8e . Labelling the 

relevant 8 graphs with variables ci, (73,(75 = ±1, the 
imaginary part of the sum of them can be written in 
the compact form 

Im V U?' r) = -rZ s (e) / d Xl dx 3 dx 5 e ~ [xi+X3+X5) 

Esm[\v\(<r 1 x 1 + a 3 x 3 + a 5 x 5 )/2] 
CTlf7 5 ZIZTHI — T^TZ. — I woi ■ ( 4 - 7 ) 



siiih[u(<7ixi + a 3 x 3 + 05X5)/ 2] 



Here, the factor a\a^ gives the signs of the respective 
graphs. 

To evaluate the remaining odd time integrations, it is 
useful to insert the identity in form of 



1 = 



ds 



dk 

e ifc[s-(a-ixi+o- 3 X3+o- 5 X5)] 

2tt 



r 



(2,6) 

















FIG. 7. The 12 graphs obtained from the 2 di- 
agrams in order Uz shown in Fig. ^. The graphs 
(1,1), (1,3), (1,5), (2,1), (2, 3), and (2,5) have positive sign, 
the remaining ones are negative. 

From Fig. [7] it is obvious that the graphs 

(1,1), (1,2), (2,1), and (2,2) fall into the class / 3 (1) , whereas 

(2) 

all other graphs belong to the class f 3 ■ To identify the 

i 3 \ 

classes giving leading logs, we put N A ' = 1 for the mo- 
ment. Then, in performing the even time integrations, 
we have to look for terms g ivin g squares of logarithms, 
v\. Upon using Eq. (4.1), it is easily shown that 



In 2 |u 



the first class /g gives subleading terms ~ In \u+v\ only, 
and solely the eight diagrams in the second class produce 
leading logs. The even time integrations are tackled by 
switching to the variables 

x = x 2 + Xi , p= (x 2 - Xi)/2 , 

which allows for exact evaluation of the p integral. As 
the logarithms come from the region x\ , x 3 , x§ x, the 
x-integral takes the same form for the eight remaining 
graphs. Upo n perfo rmin g a substitution x + q — > x as 
given in (4.2) with (4.3), where now q is a, b or c, the 
logarithms of \u + v\ are in the integral 



L 3 (e) 



dx 



In a; 

t,1-8c 



(4.6) 



-\n 2 \u + v\ (l + 0(eln\u + v\)} , 



Performing the summation over the a's and the integra- 
tion over xi,x 3 ,xs, one finds a factor — 8fc 2 /[l + fc 2 ] 3 , 
and we are left with the h- and s-integrals. These can 
be carried out analytically, and we obtain the leading-log 
terms of U 3 in the form 



Im U 3 = -4e 2 -L 3 (e) Im 

7T 



2u 



2u 



(4.9) 



2u 



This gives together with U\ and U2 all leading-log con- 
tributions to the mobility up to order e 2 . 



C. Contribution U4 

To understand the general structure of all higher-order 
terms, it is very illuminating to study the next contribu- 
tion. Although this is quite tedious, it will turn out to 
guide us towards an exact summation of all diagrams rel- 
evant for the asymptotic low-energy behavior. The five 
diagrams contributing to U4 are shown in Fig. ^. Each 
of them can be decomposed into 4! = 24 graphs, yielding 
120 graphs in total. Diagrams (I) to (IV) are of type (A), 
while diagram (V) is of type (B). Since the calculation of 
U4 is some what l engt hy bu t proceeds along the same lines 
as in Sees. IV A| and IV B, the relevant technical details 
are given in the Appendix. In this section, we only dis- 
cuss some conclusions obtained from this c om putati on. 

The main findings are exhibited in Eqs. ( [A3| ) and (A4) 
of the Appendix which show that there is a crucial dif- 
ference between diagrams of type (A) and type (B). In 
type-(A) diagrams, one has a product of only two cr's in 
the integrand, while diagrams of type (B) contain at least 
four a factors. This characteristic behavior directs us to- 
wards an important conclusion. Upon pulling out a factor 



L4(0)|w|/u 4 in Eq. (A4), we may expand the remaining 



expression in powers of u and \v\ . Now, because of the 







cr-factors, terms up to order u 2 and \v\ 2 cancel out in the 
sum over {a}, and nonvanishing contributions from the 
odd time integrations are at least of order u 4 , \v\ 4 , u 2 \v \ 2 . 
Since the asymptotic low-energy behavior in Eq. (3J) is 
governed by a,j ~ u 2 , \ v\ 2 , it is obvious that the asymp- 
totic properties are entirely unaffected by diagram (V) . It 
is immediately clear that also all higher-order diagrams 
of type (B) do not contribute to coefficients a,j ~ u 2 , \v\ 2 . 
From this we conclude that, while the zero-energy value 
of the conductance G = is fully determined by the dia- 
gram shown in Fig. |^, the asymptotic corrections in u and 
\v\ are completely governed by diagrams of type (A). Fol- 
lowing this observation, one can make further progress. 



This will be discussed in Section V A 



Let us finally collect together all contributions in order 
e 3 due to Im U4 for arbitrary pow ers in u 2 and \v\ 2 . U sing 



again an identity of the form ( |4.8| ) , the integrals ( A3 ) and 
(A4) given in the Appendix can be evaluated exactly. In 
the end, we find 



ImC/ 4 = 64e 3 -L 4 (0) Im 



2 + u + v 
2u 



(4.10) 



_ ll/ /2 + u + v\ n _ 2 ,„ /2 + u + v 



+ u - 3 ?/> (3) ( 



2u 

2 + u + v 
2u 



In view of the increasing complexity it appears to be 
prohibitive to go beyond order e 3 in all powers of u 2 , \ v\ 2 . 
However, as the expansion is already well advanced, in 
the next section the findings obtained so far will be taken 
up to conclude on exact results. 



V. RESULTS 

A. Asymptotic low-energy corrections 

Our analysis in the previous section sug gests that one 
can sum the whole leading-log series ( |3.l| ) for the con- 
tributions dj ~ u 2 , |v| 2 , whe re scaled temperature u and 



bias v are given in Eq. ( |3.4| ). These contributions deter- 
mine the asymptotic low-energy properties of the trans- 
port process. The leading-log summation can indeed be 
performed in all orders by considering only the asymp- 
totic regime. As shown in the previous section, the only 
diagrams contributing to this regime are of type (A) . The 

sum of all of them can be expressed as ( Y + +), 

where Y stands for the grand canonical gas of finite- 
length dipoles of the two types (+— ) and ( — h). (So 
Y is b asica lly a sum over all arrangements X mentioned 
in Sec. III.) In the remainder of this subsection, we will 



ignore all other diagrams .EH 

The contribution Un to U n due to diagrams of type 
(A) can now be evaluated in every order n. The even 



time integrals result in a function L n (e) containing the 
logarithm, 



L n (e) = 



1 



In" 



(n-2)! J 1 dx x 1 -^ C 



^_ -\u+v\x 



(5.1) 



Here, we have already integrated out all even time in- 
tervals except the totally symmetric combination x = 
X2 + ■ ■ -+X2n~2, and we have included the interaction fac- 



tor (3.12). Taking into account the remaining odd time 
integrals and summing over all possible arrangements of 
the intermediate dipoles, we find 



DTi- 1 

Im £ [ /0» = (-l ) ™-i i _ T i n (e) ,5.2. 

je(A),r 



: / dx\ ■ ■ ■ dx2n-i e 
10 



E 



Cl02n-1 



(xiH hX2n-l) 

sin[H(cTi2!i H h CT2n-1^2n-l)/2] 



sin 



h[«(<7iXi H h (J 2 r l -1^2n-l)/2] 



Upon expanding the integrand, we find to lowest order in 
u and 1 17 J the nth contribution to the conductance (n > 1) 



~~Gr7 



_ 2e \u + V \ 2 (-Wer : 2 



12 



(n-2)! 



In"' 2 x 



\u+v\x 



Finally, summing over all n gives to lowest integer orders 
in u and \v\ the expression 



G 



1 



Go 12 



\u + v\ 2 1 



dx- 



\u+v\x 



„l+8£ 



(5.3) 



An important feature should be noted. If one stops the 
expansion at the term n = 2, an expression of the form 
( |5.3D would emerge with the factor x 8e in the denomina- 
tor of the integrand replaced by x~ Se . The whole non- 
perturbative effect of inserting a grand canonical gas of 

finite-length dipoles between the charge pairs ( ) and 

(-) — h) simply results in exchanging the factor x~ Se in the 
denominator of the integrand for the factor x 8e , as it is 
already done in the integrand of ( |5.3|) . Clearly, this sub- 
stitution leads to a crucial change in the dependence on 
u and \v\. 

Upon performing the integral, we finally get the exact 
asymptotic behavior 



G(V,T)/G = (1/3) [(eV/2k B T K ) 2 + (nT/T K ) 



1 l+4e 



(5.4) 

This result is valid under conditions ( |2.3| ). For the lin- 
ear conductance, the T 2 / 9 ~ 2 behavior found by KF is in- 
deed recovered. More generally, the low-temperature cor- 
rections to the zero-temperature nonlinear conductance 
will follow this law as long as eV -C k-gT. However, if 
eV 3> k^T, one finds a ^-independent universal T 2 law. 
Clearly, for any finite value of the external voltage, the 
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ultimate lo w-tem perature behavior should then be T 2 . 
The result (5.4) describes a smooth turnover between 
these two power laws and demonstrates unambiguously 
that the limits V — * and T — ► do not commute. 

The nonanalytical form of the conductance is due to 
the critical nature of the (T — V — 0) model. Inclu- 
sion of a finite voltage breaks scale invariance, and if one 
considers the temperature dependence under V ^ 0, ana- 
lytical behavior will be observed. The physical reason for 
T 2 corrections is the low-frequency thermal noise of the 
plasmon modes in the leads.E3 Similar low-temperature 
T 2 features can be found in a numbauof different sys- 
tems where Ohmic damping is present.EJ We believe that 
these universal T 2 features can indeed be observed in the 
g ji v = 1/3 FQH transport experiments of Milliken et 
al.a Experiments have so far been carried out at small 
source-drain voltage eVjk^T « 0.2, and hence the KF 
T 4 law has been found. It would be interesting to ex- 
plore the regime eV ^> k^T, where one should be able 
to see the T 2 law. Sample heating is not expected to 
cause major problems since G is still quite small. Finally 
we mention a completely different experimental possi- 
bility for detection of this crossover from T 2 ^ 9 ^ 2 to T 2 
behaviors, namelw the novel quantum wire experiments 
by Tarucha et alB 



B. Conjectured exact result 

Our explicit calculation for U n up to order n = 4 in 
Sec. I V as well as the exact summation of the leading-log 



series ( |3.l| ) to lowest order u , \v\ allow us to conjecture 



on the asymptotically exact result. All these results are 
consistent with the following expression for the mobility 



G(V,T) 
Go 



1 



27r 

eV 



Im ip 



7r + ieV/2 
2vrfc B r 



(5.5) 



which is just the exact solution (3.3) for g = 1/2 with a 
voltage- and temperature- dependent renormalization 



7 -> 7r = 7 ((e^/27) 2 + (7rfc B T/7) 2 )' 



(5.6) 



We remind the reader that 7/ fee is just the Kondo tem- 
perature. Expanding Eq. (5.5) in powers of e, one finds 
exactly all contributions up to order e 3 — but in all 
orders of u 2 , \v\ 2 — as d erived in Sec. IV. Similarly, 
when exp andi ng Eq. (5.5) up to order u 2 , |u| 2 , we re- 
cover Eq. (5.4) in all orders of e. Furthermore, the T = 



limit of (5.5) -coincides with very recent calculations by 
Fendley et al£5 

It is tempting, in the light of these findings, to assert 
that the result (5.5) is the exact solution under condi- 
tions (2.3). This statement also receives support from the 
corresjpondence of our model with an interacting fermion 
model.E j To make the correspondence, we rewrite the so- 
lution (5.5) in terms of a suitable momentum integral 
representation. Upon using the representation (x and y 
real) 



Im?/;(l/2 + x + iy) 



ds 



1 -|- s 2 \ e 2it(xs-y) _|_ I 



1 



e 2-K(xs+y) _|_ I 



and putting s = (v F p/j) 1+4e , Eq. (5^) takes the form 
G 



r°° 

r =!--# / dpn(p)\S + -(p)\ 2 
Go eV J 

x (/[£(p) - eV/2] - f[E(p) + eV/2]) , 
where f(E) is the Fermi function, and where 



(5.7) 



i^wr-fi + rae) 



2+8e N 



(5.8) 



n(p) 



(v F p) 2 



2c 



(5.9) 
(5.10) 



XeV/2) 2 + (irk B T) 2 
E(p) = v F pn(p) 

In the work by Fendley et al.f3 a similar formula is 
derived for g — 1/3 using the thermodynamic Bethe 
ansatz and a Bol tzmann rate expression. The function 

1 5_| (p)\ 2 in ( |5~§| ) agrees with the form of the squared 

modulus of the exact 5 matrix of the quasiparticles for 
impurity scatteringila on putting g = 1/2 — e. The func- 
tion n(p) represents the density of states of these quasi- 
particles at voltage V and temperature T in the regime 
( |2.3| ). Finally, the quantity E(p) is a pseudoenergy used 
to parametrize the occupation probability of the quasi- 
particles as f[E{p)\. 

In conclusion, the results of our dynamical method are 
in full agreement with what one would expect from com- 
bining the thermodynamic Bethe ansatz with a Boltz- 
mann rate expression. Moreover, we obtained the den- 
sity of states and the pseudoenergy of the quasiparticles 
in explicit form for g = 1/2 — e. 



VI. CONCLUDING REMARKS 

We have studied transport properties of a Luttinger 
liquid through a barrier at low temperature and for small 
voltage. This nonperturbative problem has been ap- 
proached by a leading-log summation technique which is 
valid near the value g = 1/2 of the interaction constant. 
Putting g = 1/2 — e, the leading logs have been summed 
up in all orders e™ for the asymptotic low-energy regime. 
Furthermore, the conductance has been evaluated up to 
order e 3 for all powers of V 2 ,T 2 . Both findings allow to 
conclude on the exact result for the low-temperature con- 
ductance near g — 1/2. The asymptotic low-temperature 
corrections exhibit a smooth turnover from the T 2 / g ~ 2 
law to a universal T 2 behavior as the voltage is increased. 
In fact, for any finite voltage, one should find the T 2 
law for kftT <C eV. We have also provided a simple 
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physical argument to explain the T 2 enhancement at fi- 
nite voltage. These T 2 corrections can be found even 
in the presence of electron-electron backscattering. All 
these Luttinger liquid fingerprints should be observable 
with the fractional quantum Hall devices of Milliken et 
al.U We predict that application of a source-drain volt- 
age drop eV > feeT will change the KF T 4 law into a 
T 2 behavior. 

To conclude, we mention that our method provides the 
unique possibility of performing exact noise calculations. 
It seems very difficult or even impossible to obtain ac 
noise properties with any other technique available at 
the moment. 
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This appendix contains the technical details of the 
computation of t he con tribution [/ 4 to the conductance 
discussed in Sec. IV C| . The 96 graphs related to the 
diagrams (I) to (IV) shown in Fig. || fall into four 
classes regarding their dependence on the long even times 
X2 , X4 , xq . Denoting combinations of the short odd times 
as a, b, c, d, these four classes are given by 

fP = (sinh[u(:E2 + a)/2] sinh[u(a;4 + b)/2] 



x sinh[it(:Eg + c)/2] sinh[u(iE2 + x 4 + xq + d)/2] 
/i 2) = (sinh[w(a; 2 + a)/2] smh[u(x 6 + b)/2] 

x sinh[w(.x 2 + x± + c)/2] sinh[u(.x 4 + x 6 +d)/2] 

ff ] = fsinh[u(x 4 + o)/2] sinhKar 4 + b)/2] 

x sinh[u(.x 2 + x 4 + x 6 + c)/2] 

x sinh[u(x2 + x 4 + x 6 + d)/2]^ 
ff ] = (sinh[u(x 4 + a)/2] sinh[w(iE2 + x± + b)/2] 

x sinh[?i(x4 + xq + c)/2] sinh[w(x2 + X4 + xq + d)/2]j 



Putting = 1 for the moment in analogy to Sec. |IV B| , 
the triple integral 



dx 2 dx 4 dx 6 e -v(x2+x i+ xe) f(r) 



(Al) 
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must give the logarithms, and we are interested in terms 
~ In 3 |it + v\. It is straightforward to show that the class 

(3) 

fl does not give leading logs. The other three classes, 
however, give leading logarithmic contributions. 

Switching to relative times and to the variable x — 
X2 + X4 + the former are easily integrated out. Upon 
performing the shift x+q — ► x, where q is a, b, c, or d, and 



restoring 



(A) 



the remaining integral resulting 



from Eq. (Al) is 

1 r 



dx 



In 2 



n l-8€ 



-|fz-t-i>|a 



(A2) 



Having evaluated the even time integrals in this way, the 
remaining odd time integrals of the relevant graphs take 
forms such that the sum of them indeed gives an expres- 
sion fitting to the previous forms (4.4) and fl4.7| ), 



Im 



E 

{X1+X3+X5+X7) 



128 



{ ff =±i} 



dxidx^dx^dxj 



(A3) 



sia|H(aiXi + 03X3 + 05X5 + a 7 x 7 )/2] 
sinh[u(crixi -I- 0-3X3 + 0-5X5 + o 7 x 7 )/2] 

This contains all contributions of graphs due to the type- 
(A) diagrams (I) to (IV). Using the same procedure as 
before for U3, this expression can be evaluated without 
any further approximation. Before stating the result, we 
discuss the remaining diagram (V) in Fig. ||. 

The 24 graphs related to diagram (V) fall into three 

classes. One is just f± from above and does not con- 
tribute to leading logarithms. Another one is and 
does contribute. The third and last type is of the form 

^sinh[u(x2 + X4 + a)/2] sinh[ii(x2 + X4 + b)/2] 

x sinh[u(x4 + xq + c)/2] sinh[u(x4 + xq + d)/2]j 

and does not contribute to leading logs. In view of the 
before-mentioned complications of the even time integra- 
tions due to the residual interaction factor , we re- 
strict our attention to the order e 3 in the diagram (V), 

i.e. we put N± = 1 and do not push the explicit cal- 
culation beyond this order. Carrying out the even time 
integrations as before, we then find the factor £4(0). 

The remaining odd time integrals of the relevant 
graphs can be written in such a form that the contri- 
bution of diagram (V) emerges as 

64 /*°° 
U^' r ^ = —L±(0) / dx\dx 3 dx^dxj 



Im 



x e 



je(B),r 

{cr } =±l} 

sin[|u|(oiXi + 0-3X3 + 0-5X5 + 0-7X7V2] 
sinh [u(ctiXi + 0-3X3 + o- 5 x 5 + o 7 x 7 )/2] 



(A4) 



The crucial difference to the expression ( A3 ) for diagrams 
(I) to (IV) is the sign factor aia^a^a^ instead of o^o^. 
These distinguishing features are easily extended to 
higher-order diagrams n > 4 as follows. Every odd time 
interval T2j—i that is fenced by charges of the same sign 
gives a factor o^j-i in the integrands of the respective 

functions Im J2je(A), r U n and Im J2je(B),r u n\ see 
Eqs. ( |A3| ) and (A4). Diagrams of type (A) involve a 
product of only two a's in the integrand, while diagrams 
of type (B) imply a product of at least four cr's. The final 
result for U4 and the consequences of this observation are 



discussed in detail in Sec. IV C 
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